o 
o 

00 



On the possibility of an astronomical detection of chromaticity effects in 

microlensing by wormhole-like objects 

Diego F. Torres 1 '* Ernesto F. Eiroa 2 '} and Gustavo E. Romero 3 '''' 
1 Physics Department, Princeton University, NJ 08544, USA 
2 Institute de Astronomi'a y Ffsica del Espacio, C.C. 67, Sue. 28, 1428, Buenos Aires, Argentina 
3 Instituto Argentino de Radioastronomia, C.C. 5, 1894 Villa Elisa, Buenos Aires, Argentina 

February 7, 2008 



Abstract 

We study the colour changes induced by blending in a wormhole-like microlensing scenario with 
extended sources. The results are compared with those obtained for limb darkening. We assess the 
' possibility of an actual detection of the colour curve using the difference image analysis method. 

On t 

y—i 1 Introduction 

o : 

q This paper continues our study on the possible observational effects that struts of negative masses would 
O^roduce if they are isolated in space Q. Since wormhole structures require the violation of some of the 
5— i most sensitive energy conditions at the wormhole throat, wormholes are natural candidates -if they exist 
?Pjit all- for stellar size negative mass objects. Different wormhole solutions have been presented in the 
. £h lit era t ure after the leading work of Morris and Thorne [Q] (see for example Refs. Q). Many of these 
^ solutions actually present a negative energy density and open the possibility of having a total negative 
^ mass. However, only a few works deal with the problem of developing observational tests for the existence 
of wormhole-like objects. Our aim in the present series of papers is to turn the speculation on macroscopic 
amounts of negative masses into an experimental question, one whose answer could be reached by current 
astrophysical observations. 

In a recent paper Q, we have studied the gravitational microlensing scenario that a negative mass 
point lens would produce over an extended source. This allowed us to present more realistic light curves 
for wormhole microlensing events than those obtained earlier by Cramer et al. |J. Using the formal- 
ism introduced in Refs. ||, we computed the effects of a finite source extent on the spectral features 
of microlensing. We showed that limb darkening of the intensity distribution on a stellar source induces 
specific chromaticity effects that are very different from what is expected in the positive mass lens case. 
The possibility of using multi-colour optical observations to search for galactic or inter-galactic natural 
wormhole-like objects was then foreseen. 

Detection of the extended source effects from colour measurements, instead of single band photometry, 
is interesting because of two facts (see Ref. for further discussion). Firstly, by detecting the colour 
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curves the extended nature of the source is revealed: if the source approaches very close to the lens caustics 
but do not cross them, the induced amplification can always be mimicked by changes in the lensing param- 
eters of a point-like object. By contrast, the colour curves can not be mimicked by any such changes: a 
point source lensing event should always be achromatic. Secondly, the colour curve allows one to measure 
the lens proper motion quite easily, without the need of fitting the entire light curve. 

However, measurements of the colour curve can actually be hampered by light blending caused from 
nearby and background sources, which also causes chromaticity effects. Han et al. [|j have demonstrated 
that even for a small fraction (less than 2%) of blended light, the colour changes caused by blending can 
be equivalent in magnitude to those caused by limb-darkening. Therefore, in order to get predictions for 
a colour curve, it is essential to take blending into account, and to remove, somehow, its effects. 

In the present letter we shall analyze the chromaticity effects produced, in the case of a wormhole-like 
microlensing event, by blending of other stars. In addition, we shall estimate the likelihood of carrying 
out an actual observation of the colour curves using the difference image analysis method within current 
technological capabilities. 



2 A brief summary of microlensing formulae 

The amplification produced by gravitational lensing of a point source is given by || 

B \/Bo 2 ± 4 

where the plus sign corresponds to positive mass and the minus sign to negative mass lensing, and Bq = 
bo I Re is the lens-source separation in units of the Einstein radius Re, 



As usual, D os is the observer-source distance, D i is the observer-lens distance, Di s is the lens-source 
distance, and M the mass of the gravitational lens. For an extended circular source, instead, the light 
curve is given by j| 

A = Io n Jo* gfo <p)Ap(r, <p)rdrd<p . g . 
So" Jo* T ( r ' <f)rdrdip 

Here, (r, ip) are polar coordinates in a reference frame placed in the center of the star, r* is the radius 
of the source, and T(r, ip) is its surface intensity distribution. For a radially symmetric distribution, the 
previous expression transforms into (defining the dimensionless radius R = v/Re) 



J 2lT J ** X(R)Ap(R, <p)R dRd<p 
2nJ R * Z{R)RdR 



a= j » r r , (4) 



where R* = t*/Re is the dimensionless radius of the star. If the lens is moving with constant velocity v, 
the lens-source separation (in units of the Einstein radius) is given by 

B(T) = = J {B Q + R sin pf + (-T + R cos pf, (5) 
Re 

where T = vt/RE (see Ref. Q for a helpful plot and further details). Replacing Bo in Eq. (|l]) by its time- 
dependent partner, B(T), and using this expression in Eq. (0), we arrive, for a given intensity distribution 
1(R), at the light curves produced by lensing in an extended source case. 
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3 Blending 



The obscuration of the intensity profile of a star towards its border is known as limb darkening. An 
extended source microlensing event become chromatic as a consequence of this effect, see for instance Ref. 
||. The colour change caused by limb darkening of the source star can be computed using ||] 

A(m„ 2 -m„J = -2.51og(j^), (6) 

where A Ul and A U2 are the amplifications in two different wavelength bands, z^i and vi- For the intensity 
profile we shall take, in terms of the radius R, and as in Refs. j|, [|, 

UR) = i - C (i - /i-(f) 2 ) , (7) 

with the limb-darkening coefficients C Ul = 0.503 , C U2 = 1.050 corresponding to the I and U bands of a 
K-giant with T eS = 4750K. 



But the light curve of a microlensing event can also be chromatic by another effect: blending. Basically, 
the light flux of a source star can be affected by blended light of other unresolved stars, having themselves 
different colours, what results in a change of the colour curve. If we consider both effects at the same time, 
limb darkening and blending, the so generated colour curve was recently computed by Han et al. [R] to be 



l>2 



m 



in, 



-2.5 log 



Av 2 + fl 



I '2 



l + fu 



(8) 



where f Vi are the fractions of the blended light in the individual wavelength bands. These fractions depend 
on the specific situation and will be different for different backgrounds. The colour changes due only to 
the blending effect will be the difference between the colour curve of the limb darkening event affected by 
blending (Eq. |||) and the colour curve for limb darkening alone (given in Eq. (||) and Ref. @). To ease 
the comparison with the standard (i.e. possitive mass) case, we shall adopt the blending coefficients as 
fu = 1.39% and // = 1.04%, and a source star with radius R* = 0.1 0. Our new results, including the 
effect of blending, and for different impact parameters k = 6o/ r * = Bq/R*, are shown in Figures 1 and 2. 



The colour curves without blending present an umbra region in the negative lensing case (no light 
reaches the observer) when the impact parameter is small (k < 20) Considering blending, instead, we 
now discover that this umbra is no longer present, but rather that there is a 'plateau' (A(U — I) ~ —0.31) 
in the colour curves, produced only by the blended light. This plateau is not directly shown in the figures 
in order to show the detail in the upper portion of the colour curve. This new effect has important impli- 
cations in the full colour curves, as Figure 2 shows. 

In the case of an ordinary lens, the colour curves affected by blending are very similar to the photomet- 
ric ones, see 0] for a comparison. We see that as it gets closer to the star, the color of the observed source 
becomes redder due to the differential amplification of the coldest regions. When the lens transits towards 
the star interior, the hot center starts to dominate the amplification, producing a dramatic change in the 
slope. 

For the limb darkening colour curve (dash curve in Figure 2), the spectral changes start long before 
than in the standard situation. Initially, the source also becomes redder and then experiences a switch 
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Figure 1: Colour changes produced by blending in microlensing events. The number of blended stars 
is shown in each curve. Upper left: Positive lensing, impact parameter k = bo/r* = 0.5. Upper right: 
Negative lensing, same impact parameter. Bottom Left: Positive lensing, k = 2. Bottom right: Negative 
lensing, k = 20. The source star has dimensionless radius R* = 0.1. 
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Figure 2: Full and partial contributions to chromaticity effects in microlensing. The solid line is the full 
colour curve of an event affected by blending (two blended stars) , the dash curve represents the colour curve 
if the same event is not affected by blending, whereas the dash-dot curve gives only the colour changes 
due to the blending effect. Upper Left: Positive lensing, impact parameter k = 6o/ r * = 0.5. Upper right: 
Negative lensing, same impact parameter. The umbra is no longer present, but rather there is a plateau 
(A(C7 — I) ~ —0.31), that dominates the umbra region in the colour curves, produced only by the blended 
light. This plateau is not shown in the figure in order to show the detail in the upper portions of the colour 
curve. Bottom Left: Positive lensing, k = 2. Bottom right: Negative lensing, k = 20. The source star has 
dimensionless radius R* = 0.1. 
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when shorter wavelengths begin to dominate. Contrary to what happens with positive masses, the spectral 
trend changes again, with the source appearing colder and colder until it vanishes in the umbra during 
the transit. When the source is seen again, the inverse behaviour is observed. If we now take into account 
the blending effect as well, the existence of the previously mentioned plateau, instead of the umbra region, 
make the colour curve to change its trend again, towards the blue region. During the transit, it is the 
blended contribution the one that dominates the colour curve. It makes sense: blending fluxes come from 
stars whose light is not deflected by the wormhole-like object, and so the typical umbra effect is absent. 
Blending, then, and contrary to the positive mass case (where the pattern of the colour curve is maintained 
with only slight changes in the actual values for A(Z7 — /)), noticeably affects the form of the colour curve 
in a wormhole-like microlensing event. 

The difference between the negative and the positive colour curves (that we show for comparison in 
the same set of figures) continues to be very clear, and hence, these combined effects allow to distinguish 
between the different kind of lenses. We shall now focus on demonstrating that the colour curve can 
actually be observed with current technology in typical cases. 



4 The DIA colour curve 



The difference image analysis (DIA) is a method to measure blending-free light colour variations by sub- 
tracting an observed image from a convolved and normalized reference one. The flux would then be, within 
DIA, 

F„ = F VtOhs -F v , ie{ = {A u -l)F Vt0 , (9) 

where i^ i0 bs = A V F V $ + B v and -F^ref = F u> o + B u stand for the source star fluxes measured from the 
images obtained during the progress of the microlensing event, and from the reference (unlensed) image, 
respectively. B v is the blended flux. Then, the DIA colour curve is given by jjj] 



m 



vi J DIA 



-2.5 log 



A U2 (t) - 1\ fA U2 (t Tei ) - 1 



A ui (t) - lj \A Vl (t ief ) - 1 



(10) 



The advantage of measuring this curve, instead of that given by Eq. (||), is that it does not depend on the 
blending parameters f v j (equivalently, B v ^). We shall choose i re f from the condition 



"2 



m, 



viJDIA 



0. 



(11) 



when the reference star suffers no amplification. Basically, t re { — > oo. Again, we shall fix our attention to 
the U and I bands of a K-giant source star with dimensionless radius -R* =0.1 and T e g- = 4750-fT. The 
results for positive and negative lensing with different impact parameters k = &o/ r * = Bq/R^ are shown 
in Figure 3. 



Even when the DIA colour curve can have a different form when compared with that produced only 
by limb-darkening, they both depend on the same parameters, A Vl and A V2 . Hence, the same information 
can be extracted from both curves, but with significantly reduced uncertainties in the DIA case, because 
of the absence of blending. 

It is interesting to directly compare, then, the DIA colour curve just presented with the limb-darkening 
photometric curve presented in Figure 3b of Ref. or here in the right panels of Figure ^, dash lines. 
The analytical difference between both colour curves reduces itself to the replacement 

A A - 1 1 ' 
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Figure 3: DIA colour curves. Left: Positive lensing, impact parameter k = &o/ r * = 0.5 (solid curve) and 
k = 2 (dash curve). Right: Negative lensing, k = 0.5 (solid curve) and k = 20 (dash curve). The source 
star has dimensionless radius R* = 0.1. 




Figure 4: Left: Example of the partial evolution of the ratio (A U2 — l)/(A Ul — 1) in the case of negative 
mass lensing. Right: DIA colour curve in the same temporal interval. Impact parameter is k = 0.5, and 
R* = 0.1, for both figures. 



within the logarithm function used in the magnitude definition. This apparently simple change has, how- 
ever, large implications for the negative mass colour curve when k < 20 (6o < 2Re)- Particularly, when 
either A Ul or A U2 are less than 1, but not both, the ratio (A V2 — 1)/(A V1 — 1) is less than zero, yielding a 
not defined colour change. This happens just before the umbra, when large variations in the amplification 
suddenly occur at slightly different times for different frequencies, this being the reason of the apparent 
extra cusp in the DIA colour curve. We show the behaviour of the ratio (A U2 — l)/(A Ul — 1) for our two 
particular frequencies in Figure ^. 

Interestingly, the positive mass DIA curve is completely similar to the photometric one, since there is 
no time at which A Vl — 1 and A U2 — 1 have a different sign. 

The behavior of the negative DIA colour curve deserves further study. In order to explore exactly the 
form of the curve that could actually be measured, we would need to implement a numerical code with a 
given binning in time (corresponding to a given integration time of a telescope). If one of the cusps in the 
colour curve is produced only by a single point, we might lose it in the binning process, but we shall shed 
some light on the behavior that could actually be observed. We have then adapted the numerical code 
used in Refs. [f?J to the case of negative mass lenses. Figure [| shows two particular examples obtained 
with this code. These curves show the qualitative expected behavior in its full extent. 
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Figure 5: Negative mass DIA colour curves obtained with a numerical code, by binning the time interval 
for an impact parameter k = 10 (left) and 20 (right). The left hand side presents a break on the x-axis 
of the colour curve, where the umbra is located, in order to provide extra details. The source star has 
dimensionless radius R* = 0.1. 
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5 Measuring the DIA colour curve 



Although the colour changes are usually small, they can be measured within current technological limita- 
tions. Following Han et al. 0, we write the uncertainty in the determined source star flux as related to 
the signal-to-noise ratio by 

SF vja /F vja 1.09 
= ~0llnl0~ = W (13) 
Then, the uncertainty in the measured colour is related as well to S/N by 



1 54 

~ w . (14) 



DIA 

tional to the source flux variation, 



If S/N = 10, 5 A {m U2 — J ~ 0.15. The signal measured from the substracted image is propor- 



Soc(A v - l)F ,^ exp , (15) 



where i e xp is the exposure time. The noise comes from the lensed source as well as from the the blended 
background stars 0, 

N oc [A V F 0>V + B] 1 / 2 ^, (16) 

where B represents the average total flux of unresolved stars within a seeing disc of radius Seeing- Then, 
the signal to noise ratio is given by 

s/Ar =^- 1)F -(xra) I/2 - (17> 

Since we want to compare our error estimates with those corresponding to a positive case, we shall assume 
mutatis mutandis all parameters used in the discussion of the latter in Section 5 of Ref . . 

Let us first take the source size as 0.07 Einstein radii, and the Einstein time scale as 67.5/2 days |7|]. The 
lensed source is a K-star with 1=14.05 mag. Observations are assumed to be carried with a lm-telescope 
with a CCD camera that can detect 12 photons per second for a 1=20 mag star. The exposure, i e xp! is 
considered variable so as to allow for the measured signal to be ~ 4 10 4 photons, which is in the range of 
the linear regime response in modern CCD cameras. Actually, 



4 10 4 photons 
A v Fq >v 



''exp — 4 771 ' V-*-"/ 



and so it will be different for each given magnification. The estimation of B is done by assuming that 
blended light comes from stars fainter (i.e. with greater magnitudes) than the crowding limit, set when the 
stellar number density reaches ~ 10 6 stars deg -2 . This number density corresponds to I ~ 18.1 J?]]. The 
background flux is normalized for stars in the seeing disc with seeing = 2 arcsec. In the case of a positive 
lens, the exposure time required to achieve the requested flux of 10 photons is only about some seconds, 
and this happens due to the huge magnifications that the lensing produces (up to 20 times around t = 0). 



In the negative mass lensing situation, the overall presence of the umbra dominates part of the error 
estimation as well. In particular, for magnifications less than 1, the S/N is not well defined, since it 
becomes negative. But this happens just before the umbra, for only one point in the binned plot, and do 
not affect the correct estimate of the previously rising curve (on the left of the umbra, for instance). In 
addition, there is no sense in assigning an error to an absent signal, the umbra. We find that S/N for 
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Figure 6: Left: Error estimate in the DIA colour curve for a positive lens, impact parameter k = 0.9 and 
a source star with radius R = 0.07Re- Other parameters are discussed in the text. Right: Example of the 
error estimate in the partial evolution (left of the umbra) of the colour curve for a negative lens. Lensing 
parameters are the same as in the Left panel. 



the negative case can be around 80, with exposures times slightly larger than in the positive mass case, 
of about 8-10 s. This difference is produced by generically lower values for the magnification, which is of 
order 1, instead of the range 10-20 reached in the positive mass situation. In Figure [6| we show the case of 
k = 0.9. Note that the natural scale for microlensing, the Einstein time scale, represents half the physical 
time spread in the x-axis of the left panel in Figure ^. Then, the negative mass lens has a longer time 
evolution, since the particular peak we are showing happens already in a time scale for which almost the 
complete microlensing event occurs in the positive mass case. In Figure [/] we show the cases of impact 
parameters = 10 and k = 20, for which we have previously investigated the colour curve. Interestingly, 
due to small values of the amplification for the earliest or the latest times, the error significantly increases 
in these regions. This can be noticed particularly on the right panel of Figure [?]. Overall, it is clear, then, 
that within current observational capabilities we could be able to distinguish between ordinary and exotic 
lenses, through the analysis of gravitational microlensing chromatic effects. 

At the moment, most of the microlensing experiments do not use the DIA method in their data analysis. 
However this is already beginning to change, see for instance Ref. 1|, and will become a common practice 
in the near future. If the microlensing alert systems are adapted to take into account the possible colour 
and light curves produced by negative mass lenses, we shall be in position to make extensive searches -and 
to establish bounds on the possible existence- of wormhole-like objects. 
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Figure 7: Error estimate in the DIA colour curve for negative lenses, impact parameter k = 10 (again 
showing a break, left) and 20 (right), and a source star with radius R = O.IRe- 
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6 Concluding remarks 



All theoretical constructs thought to represent features of the real world should be queried through exper- 
imental or observational tests. This process is fundamental for science. In this paper, we have expanded 
the formalism for wormhole-like gravitational microlensing of extended sources by including the analysis 
of the effects of blending. Having so constructed a complete colour curve, taking into account the effects 
of limb darkening as well, we analyzed the possibilities for an actual detection of chromaticity effects. 

Struts of negative masses, if they exist at all, will be detected through the effects they produce upon 
the light coming from distant sources. If a consistent lensing survey yields a negative result, we could then 
set empirical constraints from a statistical point of view to the amount of negative mass in the universe. 
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